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ABSTRACT 

The cluster abundance and its redshift evolution are known to be a powerful tool to constrain the amplitude of 
mass fluctuations erg and the mass density parameter 57„,o. We study the impact of the finite decay rate of cold dark 
matter particles on the cluster abundances. On the basis of spherical model in decaying cold dark matter universe, 
we calculate the mass function of clusters and compare it with observed cluster abundance. We find the decay 
of cold dark matter particles significantly changes the evolution of the cluster abundance. In particular, we point 
out that the lifetime of dark matter particles comparable to the age of the universe lowers the ratio of the local 
cluster abundance to the high-redshift cluster abundance and can account for the observed evolution of the cluster 
abundance quite well. The strong dependence of the cluster abundance on the decay rate of dark matter suggests 
that distant cluster surveys may offer clues to the nature of dark matter. 

Subject headings: cosmological parameters — cosmology: theory — dark matter — galaxies: clusters: general 
— large-scale structure of universe 



> 

o : 

(N . 
O ■ 

\o ■ 
o : 

m . 

o ■ 

Or 
i 1 

o ■ 
a: 

C/J . 



1. INTRODUCTION 

The abundance of cluster of galaxies places a strong con- 
straint on fundamental cosmological parameters such as the 
amplitude of mass fluctuations on a scale of 8/r'Mpc, a%, and 
the present mass density parameter, f2 m o (Henry & Arnaud 
1991; Bahcall & Cen 1992; White, Efstathiou, & Frenk 1993; 
Eke, Cole, & Frenk 1996; Kitayama & Suto 1997). Although 
the abundance of clusters is strong function of both fi m o and 
Us, the evolution of the cluster abundance with redshift breaks 
the degeneracy between fi m o an d eg (Eke et al. 1996; Henry 
1997, 2000; Bahcall & Fan 1998; Bahcall & Bode 2003). For 
instance, Bahcall & Bode (2003) pointed out that the relatively 
high abundance of massive cluster at z > 0.5 prefers large mass 
fluctuations erg = 0.9- 1. On the other hand, the local cluster 
abundance in both optical (Bahcall et al. 2003) and X-ray (Ikebe 
et al. 2002; Seljak 2002) wavebands suggests low mass fluctu- 
ations er s = 0.7-0.8 for fl m Q ~ 0.3. The combination of local 
and high-z cluster abundances yields relatively low matter den- 
sity in the universe, e.g., f2 m o = 0.17 ±0.05 (Bahcall & Bode 
2003). However, other observations, such as the anisotropy of 
cosmic microwave background (Spergel et al. 2003) and Type 
la supernovae (Tonry et al. 2003), prefer somewhat larger mat- 
ter density; Q m0 = 0.27 ± 0.04 from WMAP plus large scale 
structure (Spergel et al. 2003) or f2,„o = 0.28 ±0.05 from Type 
la supernovae (Tonry et al. 2003). Although all these results are 
consistent with each other within 2er level, the slight difference 
of erg or f2 m o might be interpreted that underlying models are 
wrong (Bridle et al. 2003). 

Although the above constraints are derived assuming the 
usual stable cold dark matter model, the evolution may be con- 
siderably different if cold dark matter particles are not perfectly 
stable and decay into relativistic particles with a decay rate T 
(Cen 2001). Therefore, the evolution of the cluster abundance 
may become a useful test to explore the nature of dark matter. 
The decaying cold dark matter model has been studied so far 
in several contexts, e.g., to reconcile Einstein-de Sitter universe 
with current low matter universe (Turner, Steigman, & Krauss 
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1984), to explain unexpected X-ray observations (Dilella & 
Zioutas 2003), or to reduce possible over-concentration of dark 
matter and over-abundance of substructures in cold dark mat- 
ter model (Cen 2001; Sanchez-Salcedo 2003). Moreover, cold 
dark matter with similar property, disappearing dark matter, has 
also attracted much attention in the context of the brane world 
scenario (Randall & Sundrum 1999a,b). In the brane world sce- 
nario a bulk scalar field can be trapped on the brane, which 
is identified with our universe, and become cold dark matter. 
But since the scalar field is expected to be metastable, it decays 
into continuum states in the higher dimension with some decay 
width r, which is determined by the mass of the scalar field and 
the energy scale of the extradimension (Dubovsky, Rubakov, & 
Tinyakov 2000). An observer on the brane sees this as if the 
scalar field disappears into the extradimension leaving some 
energy on the brane. This energy behaves in the same way as 
energy of relativistic particles but does not correspond to real 
particles. Thus this is called "dark radiation". Although the 
disappearing dark matter and the decaying dark matter are not 
the same, they contribute to the expansion of the universe in the 
same way. Hereafter we call both of them decaying dark mat- 
ter. Ichiki et al. (2003) studied the universe with the decaying 
dark matter as a test of brane world scenario. They found that 
decaying cold dark matter model with lifetime about the age of 
the universe improve the fit of Type la supernovae observations 
and the evolution of mass-to-light ratios of clusters. 

In this paper, we study structure formation in the decay- 
ing cold dark matter universe; we present the spherical model, 
make theoretical predictions for the cluster abundance, and 
show that the evolution of the cluster abundance is quite sensi- 
tive to the decay rate of dark matter particles. We also compare 
our predictions with observed low-z and high-z cluster abun- 
dances, and see if the decaying cold dark matter model explains 
the relatively large cluster abundance at high-z universe. 



2. THEORETICAL MODEL 
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To study the cluster abundance, we need the mass function 
of dark halos in decaying cold dark matter model. In this sec- 
tion, we construct the mass function using the spherical model 
as in Press & Schechter (1974). We also present the nonlin- 
ear overdensity A c which is important in studying the cluster 
abundance. Throughout this paper we assume a flat universe. 

2.1. Cosmology and Linear Growth Rate 

We consider the case that dark matter particles decay into rel- 
ativistic particles, which are assumed to contribute to the struc- 
ture formation only through the change of the expansion law of 
the universe by their energy behavior. Rate equations of matter 
and radiation components in the decaying dark matter model 
are p m + 3Hp m = -Yp m and p r + 4Hp r = Tp m . This results in the 
following evolution of matters and radiations: p m = p m oa~ 3 e~ r ' 
and p r = Tp m oa~ 4 j^ae^'dt, respectively. Here we assume no 
radiation component except for decay products. The Friedmann 
equation then becomes 

^ = n m0 a- 3 e- r ^ + m mQ a- 4 f ae^'^dt + X,. (1) 
#0 Jo 
This equation should be solved with the equation of the cosmo- 
logical time t = L da/ (aH(a)) simultaneously. By combining 
these two equations, we obtain a following differential equation 
H 2 /O n p- r< -'-'o'> \ t' 

^" = 4 ^2 + 4 fl A ) f ' 3 A (2) 

2 \ a z J a 

where ' = d/da. Note that the current density parameter fi m o 
should satisfy the Freidmann equation at z = 0: 



i = n mQ + m mQ [ ae- T{ '-' 0) dt + A . 
Jo 



(3) 



Therefore the current density parameter £l m o is uniquely deter- 
mined for given T and Ao, i.e., f2 m o = £l m o(T,Ao). Next con- 
sider the linear perturbation. The density perturbation is de- 
noted by S m = (p m - p m ) /p m , where p m is background matter 
density. Then perturbation equation becomes 

a 2 Sy t'm-^Gp m 5 m = 0. (4) 

We solve this differential equation with the boundary condition 

S m oc a at a <C 1 . 

2.2. Spherical Model 

In this subsection, we describe the spherical collapse model 
in decaying dark matter model and give the methodology to cal- 
culate nonlinear overdensity A c . First we consider a spherical 
overdensity with initial mass M and radius R. The equation of 
motion of such a spherical shell is given by 

R _ GM 

We define the following quantities: x = a/a t . d , y = R/R ti , 
t = H(x= lWU m (x=l)t, 7 = T/(H(x = l)V$2m(jc= 1)), and 
r] = X(x = =1). Here a Vd represents the scale factor at 

turnaround. Then equations we should solve are the two di- 
mensionless equations: 



-rt 



-H^rn mQ a- 4 [ ae- r( '-' o) dt + Hr]\o. (5) 
Jo 



dx 2 



e -7(T-T la ) 



- + 4t]x 



drY 



1 dr 



dx J x dx 1 
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dx 2 



d 2 r 
dx 2 



dy 
dx 



drY 
dx J 



where £ is the density contrast of the spherical overdensity 
region at turnaround, ( = /O sp hericai/Pm|^i- We solve these 
equation from r = to the virialization of the overdensity re- 
gion with boundary conditions t(x « 1)~ (2/3)e" 7T,! " /2 x 3 ' /2 , 
t(x = 1) = T ta , y(x = 0) = 0, y(x = 1) = 1, and dy/dx(x = 1) = 
in order to derive r ta and (. Although the conventional choice 
of the virialization epoch is the collapse time (y = 0), in this pa- 
per we define the virialization epoch when the virial theorem in 
the following holds: 



dy_ 
dr 



-7(t-t„) _ 7} ,2 x -4 



/ xe-^-^dr + rjy 2 . (8) 
Jo 



Given the virialization epoch and virial radius, the nonlinear 
overdensity can be calculated from A c (x v ; r ) = (x 3 ir /y 3 ir )C- 1 . To 
calculate the cluster abundance, we also need the extrapolation 
of the linear fluctuation S c at virialization. We use S c = 1.58, 
which is the value in Einstein-de Sitter universe for all cos- 
mological models, mainly because of the computational cost. 
This simplification does not change our results because we con- 
firmed that S c only weakly depends on the cosmological param- 
eters. 



2.3. Mass Function 

We derive the mass function of clusters using the theory of 
Press & Schechter (1974). First start from an initial density 
field S(x,Mi,Zi) smoothed over the region containing mass M,. 
If the initial density field is random Gaussian, a probability dis- 
tribution function of 5 at any point is given by 



P[5(M i , Zi )] = 



1 



(2tt) 1/2 a Mi (zd 



exp 



S 2 (M h zd 



(9) 



where cr M .(zi) = cr(RM n Zi) is the mass variance. From the spher- 
ical model, it can be interpreted that the region is already virial- 
ized at z if the linearly extrapolated density contrast <5 linear (M, ,z) 
exceeds the critical value S c . Therefore the probability that the 
region with mass M is already virialized is given by 



(10) 



where erfc(x) is the complementary error function, S c (z) = 
S c D(z = 0)/D(z), and <j Mi = o~m,(z = 0), where D(z) is linear 
growth rate calculated from equation (4). We assume the mass 
variance for the cold dark matter fluctuation spectrum with the 
primordial spectral index n = 1, and adopt a fitting formula pre- 
sented by Kitayama & Suto (1996). From equation (10), we 
finally obtain the comoving number density of halos of mass M 
at time z: 



dn 



PS 



dM 



(M,Z): 



„2n 



2 po S c (z) 
ir Mi a 2 



d&\ 



dMi 



exp 



2<. 



-\y- 2 Ce-^-iyx- 4 J\e-^ , (7) withpo = Pcnt ( z = o)O m0e -n-o). 



M,=Me r ' 

(11) 
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3. EVOLUTION OF THE CLUSTER ABUNDANCE 

In the model described in the previous section, the abundance 
and its redshift evolution of clusters are fully determined by 
three independent parameters; Ao, as, and T. In this section, we 
see the dependence of the cluster abundances on these parame- 
ters. 

Figure 1 plots the number density of clusters as a function of 
redshift. We consider the cases with the lifetime of dark mat- 
ter F" 1 = 10/r'Gyr, and also with the usual infinite lifetime of 
dark matter for reference. Following Bahcall & Bode (2003), 
we adopt the mass within a comoving radius of 1 .5/z _1 Mpc, 
Mi 5h -i Mpc , to compare our theoretical predictions with observed 
cluster abundance. For the extrapolation to 1.5/z _1 Mpc, we use 
the observed cluster profile, M(< R) ~ R 6 (Carlberg, Yee, & 
Ellingson 1997; Fischer & Tyson 1997). For the observed clus- 
ter abundance, we plot the data compiled by Bahcall & Bode 
(2003); abundance at z ~ 0.05 from Ikebe et al. (2002), at 
z ~ 0.38 from Henry (2000), at z ~ 0.5-0.65 and z ~ 0.65-0.9 
from Bahcall & Fan (1998). Obviously, inclusion of finite life- 
time greatly improves the situation: for F" 1 = 10/z _1 Gyr both 
local and high-z cluster abundance can be reproduced with the 
same <j%. For Y~ l = oo it is difficult to explain local (z ~ 0.05) 
and high-z (z ~ 0.3-0.9) cluster abundances simultaneously, al- 
though both can be reconciled if we adopt larger Ao, Ao ~ 0.8. 
This drastic change of the evolution of the cluster abundance 
implies that it can place a strong constraint on the decay rate of 
cold dark matter. 
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FIG. 1. — Evolution of the cluster abundance for Ao = 0.7. The number 
of clusters with mass M l 5^-iMpc > 8 x 1O 14 /t'M0, where M i 5A -i Mpc is the 
mass within a comoving radius of 1.5/i~'Mpc, is plotted as a function of red- 
shift. The filled circles with errorbars are observed abundance (see text for 
details). Solid lines are predicted cluster abundance with lifetime of dark mat- 
ter r~' = lO/r'Gyr, while we assume no decay for the dotted lines. We plot 
curves with erg = 0.7, 0.8, 0.9, 1.0, 1.1, and 1.2 (only solid lines are labeled). 
Since adopted mass function in this paper is based on the theoretical model of 
(Press & Schechter 1974) which is now known to be inaccurate, curves may be 
different in more accurate calculations. In particular, our mass function may 
underpredict the large-mass halos, and this may result in faster decline of the 
curves with z. 

To see how large decay rate is required to explain the ob- 
servations, in Figure 2 we show the constraints from low-z 
(Bahcall et al. 2003) and high-z (Bahcall & Bode 2003) clus- 
ter abundances assuming A = 0.7. We find that the lifetime 
of T" 1 < 10/z _1 Gyr can explain both low-z and high-z cluster 



abundances. Interestingly, this decay rate is roughly consis- 
tent with that derived from the evolution of mass-to-light ratios 
T" 1 ~ lO^Gyr (Ichiki et al. 2003). 
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FIG. 2. — Constraints from high-z (solid) and low-z (dotted) cluster abun- 
dances in the r~ l -a$ plane. The observed abundance of high-z clusters is taken 
from Bahcall & Bode (2003), and that of low-z clusters is from the SDSS clus- 
ters (Bahcall et al. 2003). Contours of 68% and 95% limits are plotted. We 
adopt A = 0.7. 

We plot the constraint in the (1 -Ao)-<T8 plane in Figure 3 
with T" 1 = 7/z _1 Gyr. We also show the constraint from the ob- 
servation of Type la supernovae complied by Tonry et al. (2003) 
assuming the same value of T" 1 . We use the sample of 172 su- 
pernovae with z > 0.01 and Ay < 0.5 mag to avoid possible 
systematics. The resulting limit on 1 - Ao is 1 - Ao = 0.22^ 
(68% confidence limit), slightly lower than the limit in the case 
of no decay (1 — Ao ~ 0.28). We find that the best-fit parameter 
set (1 - Ao, as) = (0.2, 1 .2) is consistent with the result of Type 
la supernovae as well as the evolution of the cluster abundance. 




1-X 

FIG. 3. — Constraints in the (1 - \o)-a& plane assuming r -1 = 7/r'Gyr. We 
also plot the 68% and 95% confidence limits from 172 Type la supernovae 
complied by Tonry et al. (2003). The best-fit parameter set of the combination 
of these three constraints is (1 - Ao,<r 8 ) = (0.2, 1.2). 
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4. SUMMARY 

We have presented structure formation model in decaying 
cold dark matter universe. More specifically, we have calcu- 
lated the spherical model and construct the mass function in the 
similar way as the theory of Press & Schechter (1974). Us- 
ing this mass function, we have shown that decaying cold dark 
matter model with T" 1 < 10/z -1 Gyr can explain observed cluster 
abundances at both low-z and high-z universe even better than 
the usual stable cold dark matter model. Our model may be also 
consistent with the excess of cosmic microwave background 
fluctuations on small scales; large mass fluctuations erg <~ 1.1 
are required to account for this excess if it is interpreted as the 
Sunyaev-Zel'dovich effect of the cluster at z <~ 1 (Komatsu & 
Seljak 2002; Bond et al. 2003). 

While we have concentrated on the evolution of the clus- 
ter abundance, it is interesting to ask whether the decaying 
cold dark matter model is consistent with other observations 
such as the anisotropy of cosmic microwave background and 
large-scale structure of the universe. For instance, Perci- 
val et al. (2001) constrained the cosmological parameters as 
fl m ()h = 0.20 ± 0.03 from the power spectrum of galaxy clus- 
tering. Since the power spectrum mainly measures the epoch 
of the matter-radiation equality, Q m oh should be interpreted as 
Qmoe r '°h in the decaying cold dark matter model. In the case 
of h = 0.7, T" 1 = 7/z _1 Gyr, we constrained 1 - A as 0.19 < 
1 - Aq < 0.23 with 68% confidence level. This allowed range 
corresponds to 0.23 < Q mQ e r '°h < 0.27, thus marginally consis- 
tent with the observed power spectrum. Moreover, the situation 
becomes better if the Hubble constant is smaller than 0.7. For 



the anisotropy of cosmic microwave background, it is difficult 
to infer the consequence of decaying cold dark matter model 
because many parameters degenerate. We plan to study com- 
prehensive constraints of these observations on decaying cold 
dark matter model, which will be presented elsewhere. 

We note that our modeling of mass function may be inaccu- 
rate because it has been claimed that the mass function of Press 
& Schechter (1974), particularly at high mass region, tends to 
underestimate the abundance of the massive halos compared 
with iV-body simulations (e.g., Jenkins et al. 2001). This may 
be a part of reason that the cluster abundance in our modeling 
based on the theory of Press & Schechter (1974) seems to de- 
cline slightly faster with redshift than that in Bahcall & Bode 
(2003) who used the result of TY-body simulations. Thus the 
constraints on the decay rate T (e.g., Figure 2) may be weaken 
a bit if we use more accurate mass function in decaying cold 
dark matter model. One way to refine the model is to resort 
TY-body simulations, as in the case of the standard cold dark 
matter model. Nevertheless we believe that our qualitative re- 
sult that decaying cold dark matter model tends to lower the 
ratio of the local cluster abundance to the high-redshift cluster 
abundance is still unchanged even if we adopt more accurate 
mass functions. The cluster abundance at higher redshift uni- 
verse z ~ 1-2, which will be obtained by future surveys of 
X-ray or Sunyaev-Zel'dovich effects, can tightly constrain the 
decay rate of dark matter. 

We thank Kiyotomo Ichiki for useful discussions. K.T.'s 
work is supported by Grant-in- Aid for JSPS Fellows. 
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